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Abstract 

Molecular-dynamics simulations of fracture in metallic glass-like systems are 
observed to undergo embrittlement due to a small change in interatomic po- 
tential. This change in fracture toughness, however, is not accompanied by 
a corresponding change in flow stress. Theories of brittle fracture proposed 
by Freund and Hutchinson indicate that strain rate sensitivity is the control- 
ing physical parameter in these cases. A recent theory of viscoplasticity in 
this class of solids by Falk and Langer further suggests that the change in 
strain rate sensitivity corresponds to a change in the susceptibility of local 
shear transformation zones to applied shear stresses. A simple model of these 
zones is develped in order to quantify the dependence of this sensitivity on 
the interparticle potential. 

I. INTRODUCTION 

This paper presents simulations in which a small change in interparticle potential leads 
to a qualitative change in ductility. Section |IJ describes the simulation technique and ob- 
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servations. Section [II J details a technique for calculating a quantitative measure of local 
non-amne deformation which is applicable to materials with no crystalline order, and using 
this technique pinpoints those areas of the material which are undergoing some molecular 
level rearrangement akin to a dislocation. Section [IV] relates this observed change in ductility 
to the particular change in interparticle potential. This final section discusses the simula- 
tions in terms of current phenomenological theories of brittle and ductile behavior due to 
Freund and Hutchinson0 and a theory of this mechanism of molecular level rearrangement 
in non- crystalline materials developed by Falk and Langer which we shall refer to as FL.i 
Finally, a simplified microscopic model is analyzed. This model directly relates the interpar- 
ticle potential to the parameter which controls deformation in FL and the change of fracture 
toughness observed in the simulations. The model also explains the observations of Srolovitz, 
et. al. relating these regions to "r-defects" in previous metallic glass simulations.i'§ 

The concepts of brittleness and ductility are central to any understanding of failure in 
solids. The most developed first-principles theories of ductility are rooted in the dynamics 
of dislocations in crystalline solids.i~i Although it has been conjectured that an analog to 
a crystal dislocation exists in noncrystalline solids,! it remains unclear how to make the 
direct connection to molecular level phenomena necessary for these theories to be useful in 
quantitatively understanding transitions between brittle and ductile behavior in disordered 
materials. In fact, it is not at all clear that a dislocation model of this sort is the most 
appropriate way to understand materials without regular structure although these materials 
are observed to undergo similar brittle-ductile transitions to their crystalline counterparts. 
It is conjectured here that dislocation concepts are not the most natural way to describe non- 
crystalline solids and theories of "shear transformation zones" (STZ's) as first considered in 
the literature on metallic glassesi~0 are developed as a natural way to understand some of 
the basic physics of brittle versus ductile fracture. 
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II. SIMULATIONS OF BRITTLE AND DUCTILE FRACTURE IN A 

NONCRYSTALLINE SOLID 



This section describes a series of molecular-dynamics(MD) simulations of fracture in a 
simple, two-dimensional amorphous solid. While these simulations are clear examples of 
change in ductility induced by a change in interparticle potential the important point here 
is not simply to differentiate between brittle and ductile behavior, but rather to establish 
a connection between this particular change in the potential and a change the observed 
fracture toughness. 

From a practical standpoint these studies of brittle and ductile behavior are relevant in 
the context of several different disordered and amorphous materials. The simulated system 
is most similar to metallic glasses which have been observed to undergo transitions between 
ductile and brittle behavior both as a function of temperature and due to small amounts of 
dilute crystallization produced during annealing.EUll Similar transitions are also critical to 
the processing of colloidal ceramic systems. These clay-like materials undergo brittle-ductile 
transitions due to changes in salt content, i.e. changes in interparticle interactions.^ Issues 
of brittleness and ductility are also crucial for the production of high-strength polycrystalline 
metallic alloys in which such transitions have been studied experimentally with respect to 
temperature and loading rate.0 

In modeling interparticle interactions a Lennard- Jones (L J) potential is employed. This is 
consistent with previous investigations which have been carried out in the context of metallic 
glasseslOHiil! While more sophisticated models of interactions within metals exist many, 
such as the Johnson potential for ironil or the model potential for copper and zirconium 
employted by Deng, Argon and Yipjl have a basic form similar to that of the LJ potential. 
The use of a simplified potential can be justified here since these simulations seek to elucidate 
how a well controlled change in potential affects fracture behavior rather than exploring the 
accuracy of a particular potential. 

Despite the simplicity of the LJ potential, the similarity of the results of these simulations 
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to experiments carried out in metallic glasses is striking. First, as in the simulations, both 
brittle and ductile fracture are observed in metallic glasses at low temperature depending on 
composition. Typically these glasses display a pseudo-cleavage fracture mode which involves 
significant flow at the crack tip evidenced by vein patterns.il However, compositional changes 
can lead to brittle modes of fracture in which such flow is not in evidence.il Some Pd-Si 
glasses have even been observed to display the "mirror" , "mist" , "hackle" behvior typical of 
the extremely brittle oxide glasses.il Though a transitional temperature to this brittle mode 
is seen in many metallic glasses, some have been observed to experience pseudo-cleavage at 
all observed temperatures even as low as 76K.il 

Secondly, as in the simulations, metallic glasses at low temperature are observed to ex- 
perience changes in fracture toughness independent of flow stress. Tests of the dynamic 
fracture response and the onset of flow in Pd-Si and Fe-P-C glasses show a distinct crossover 
from a thermally activated to an athermal mechanism for fracture and flow at temperatures 
of 273K and 473K respectively.0'il In the high temperature regime flow and fracture are 
observed to follow the same trend when strain rate is varied. But, in the low tempera- 
ture regime analogous to these simulations, the flow stress is observed to be independent of 
strain rate while the fracture stress is strain rate dependent.il In these simulation a change 
in potential produces a change in fracture toughness that is not accompanied by a corre- 
sponding change in the flow stress. This can be related theoretically to a change in strain 



rate sensitivity independent of the flow stress, and will be discussed in detail in Section IV 



A. Methodology 

The simulated systems consisted of 90,000 particles in two-dimensions interacting via 
a two-body potential. In order to avoid problems of local crystallization, a poly-disperse 
collection of particles was simulated. The system was composed of eight different species in 
equal proportion with radii r 1; r 2 , . . . r 8 such that 
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r a = l.lr a _ x ; Yl = ^(d /2) 2 (2.1) 

0=1 

Thus, the total volume of the collection was the same as if the particles were all of radius 
do/ 2, and the system is roughly comparable to a single component system in which the rest 
spacing between two molecules is do. All quantities will be given in terms of dimensions for 
which do is the length scale. Therefore in these units do = 1. The masses of all particles 
were taken to be m = 1 in these units. 

The inter-molecular potential was different in the two simulations. In the simulation 
which displayed more ductile behavior the potential was a standard Lennard- Jones 6-12 
potential 




where r is the interparticle distance and r a and rp are the radii of the two particles, e, 
the depth of the energetic minimum of the two particle interaction, is unity in our units. 
Note, this is not the standard form of writing the LJ potential, which is typically written 
in terms of the hard-core radius, This expression is, however, equivalent and 

will facilitate the introduction of the second potential to which it will be compared. 

The simulation which displayed more brittle behavior employed a potential which will be 
refered to as the Compressed Lennard- Jones(CLJ) potential because it is the Lennard- Jones 
potential rescaled around the center of the potential well, 

U^ J {r) = U^p (A r + (1 - A) (r Q + 77,)). (2.3) 

Note that the standard LJ potential is recovered when A = 1. Furthermore, A — * cor- 
responds to a mean field limit in which every particle interacts with every other particle 
equivalently, and A ^> 1 is the limit of solely nearest neighbor interactions. For the second 
simulation the parameter A was chosen to be 1.5. This means that width of the potential well 
was smaller by 33%, and, consequently, the effective range of interaction was also shortened 
compared to the standard LJ interaction. For the sake of comparison Fig. [I] shows both 
potentials. In both cases interactions were cut off at a range of d c ~ 2.2do- 
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All times are given in units of to — doym/e. This unit of time is approximately equivalent 
to one molecular vibrational period. 

The initial amorphous systems were created by taking 10,000-molecule systems and equi- 
librating them using a sequential MD algorithm with periodic boundary conditions, a Nose- 
Hoover thermostat&0 and Parrinello- Rahman barostat JH1 The equations of motion and 
time constants for the thermostat and barostat were the same as those in FL.I The time 
step in the simulations was taken to be O.Olto- The systems were held at low temperature 
kT = O.Ole for 5000 timesteps at zero pressure, then the pressure was raised to 10e/dl over 
the course of 1000 timesteps and lowered again over an equal period of time. Subsequently 
the samples were allowed to equilibrate at zero pressure for 1000 timesteps. This procedure 
created close-packed samples. The LJ sample was observed to have a Young's modulus of 
34 and a shear modulus of 10; the CLJ sample was observed to have a Young's Modulus 
of 39 and a shear modulus of 12. These 10,000-molecule samples were then used to create 
larger systems by replicating the small system in a 3 x 3 array. 

The larger systems were simulated via a parallel MD algorithm based on a spatial de- 
composition method.0 In order to create the initial conditions for the fracture simulations, 
the large system was equilibrated for 100 time steps while held at a very low temperature, 
kT = O.OOle. A crack was then introduced into the sample. This was accomplished by 
imposing displacements as determined by the analytical solution for a straight crack in an 
elastic medium loaded below the ideal critical stress. The faces of the crack were marked so 
that the top face would not interact with the bottom face to prevent the crack from healing. 
The outer boundaries of the system were held fixed while the simulation was again run for 
20,000 timesteps holding the temperature constant to allow the system to relax. 

In the fracture simulations no thermostat or barostat was employed. To drive the crack, 
an initial velocity gradient was imposed across the sample, and the top and bottom surfaces 
were constrained to move apart vertically such that the side closer to the crack would 
separate at a strain rate of O.OOOltg 1 an d the side farthest from the crack would not move 
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apart at all. The horizontal motion of these surfaces was unconstrained. The left and 
right surfaces were constrained not to move in the horizontal direction, though their vertical 
motion was unconstrained. 

A strain rate of O.OOOli^ 1 corresponds to a physical strain rate on the order of 10 8 s _1 . 
While this may seem high compared to typical laboratory values, the time for the stresses 
to equilibrate is equal to the time for a sound wave to traverse the sample, ~ 300to- The 
fact that this time multiplied by the strainrate is ~ 0.03 1 implies that the system was 
loaded nearly quasi-statically. That is to say that the loading rate was much slower than 
the elastic response time, although the loading may not be slow when compared to the time 
scale for plastic response. Of course, if the crack begins to propagate strain rates near the 
tip may be significantly higher. 



Figure shows the average stress measured during both simulations. In order to better 
compare the two systems, the stresses are given in units of the critical stress for initial failure 
of a perfectly brittle solid with the identical elastic properties, 



Here we assume that the system can be treated as a crack in an infinite medium. E is 
Young's Modulus; a is the initial length of the crack; and G is the energy release rate, which 
can be expressed as a surface energy and a dissipation per unit crack extension, 



For an ideally brittle solid all the elastic energy released goes into the creation of new surface, 
Gdiss = 0. 7 was measured by taking a sample of the material in MD, slicing it along an 
arbitrary plane and measuring the change in potential energy. The value of 7 is 1.04e/c?o m 
the CLJ system and 0.94e/do in the LJ system, thus a l c deal = 0.68e/dl in the ductile system 
and 0.70e/d,Q in the brittle system. 



B. Observations 





G = 27 + G diss . 



(2.5) 
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Two notable differences are observed between the simulations: (i) In the CLJ simulation, 
some modest amount of energy was dissipated and the crack began to propagate at about 
7% above the ideal brittle critical stress, but in the LJ simulation fracture did not proceed 
until the stress was 48% above this value. This means that for the CLJ case the ratio of 
energy dissipated to the energy expended creating surface is 0.14, while for the LJ case this 
ratio is 1.19. (ii) In the CLJ simulation, once the crack began to propagate, the stress in 
the system sharply dipped as the crack moved through the system at speeds reaching 30% 
of the shear wave speed. Throughout this process the crack tip remained atomically sharp. 
The process stopped short of releasing all the stress because the crack arrested. In the LJ 
case, however, the crack tip blunted significantly. In this simulation, the stress remained 
high while voids nucleated ahead of the tip. The speed of the ductile crack, while difficult 
to measure due to the mechanism of propagation, stayed well below the speed of the CLJ 
crack. 



III. QUANTIFYING LOCAL DEFORMATION 

The simulation which utilized the CLJ potential resulted in markedly more brittle be- 
havior than the simulation that utilized the LJ potential. In order to address why this par- 
ticular change in potential resulted in differing amounts of deformation near the crack tip 
the underlying mechanism of deformation must be established. Work by Argon and Spaepen 
suggests that localized deformations, or "shear transformation zones," are responsible for 
rearrangements in these amorphous materials. This section undertakes an examination of 
the microscopic nature of the plastic rearrangement in order to determine if this is indeed 
the case for this set of simulations. These microscopic observations also serve to differen- 
tiate the simulations performed here from similar investigations of fracture undertaken in 
crystals@S where plasticity is observed to result from dislocations emitted from the crack 
tip or activated in the vicinity of the crack. 



S 



A. Definition of D 2 , 



In a perfect crystal, dislocations can be readily identified by their characteristic stress 
fields or as regions of anomalously high potential energy. In glasses, however, such analyses 
are difficult due to inhomogeneities frozen into the structure. Furthermore, it is not clear 
that any analog of crystalline dislocations exists in non-crystalline solids. For these reasons, 
a different scheme must be developed to identify regions which deform in a non-affine way 
and thereby observe what sort of microscopic structures play the role of dislocations in these 
materials. 

To identify local rearrangements from a set of molecular positions and subsequent dis- 
placements the closest possible approximation to a local strain tensor is computed in the 
neighborhood of any particular molecule. The neighborhood is defined in this case by the 
interaction range, d c . The local strain is then determined by minimizing the mean square 
difference between the the actual displacements of the neighboring molecules relative to the 
central one, and the relative displacements that they would have if they were in a region of 
uniform strain Sij. That is, we define 

D 2 (t, At) = EE ['»(*) - rj(f) - E(*y + ^){<{t - At) - r>(t - At))} 2 , (3.1) 

n i j 

where the indices i and j denote spatial coordinates, and the index n runs over the molecules 
within the interaction range of the reference molecule, n = being the reference molecule. 
r l n (t) is the i'th component of the position of the n'th molecule at time t. We then find the 
Eij which minimizes D 2 by calculating: 

X i3 = E K{t) - ri(t)) x (rl(t - At) - r> (t - At)), (3.2) 

n 

Ya = E K{t - At) - r'(f - At)) x K{t - At) - r> (t - At)), (3.3) 

n 

e ij = ZX ik Yj c 1 -8 ij . (3.4) 
k 

The minimum value of D 2 (t, At) is then the local deviation from afflne deformation during 
the time interval [t — At, t]. This quantity shall be refered to as -D^m- 
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B. Molecular Level Observations 

D 2 min serves as a diagnostic for identifying where local rearrangements have taken place. 
The right-hand frames in Figures [| and ^ are shaded by the value of over the interval 
from t = to the current time. It is immediately apparent that much more non-affine 
rearrangement takes place in the LJ simulation than in the CLJ simulation. In addition, 
there seem to develop preferred directions along which deformation takes place. These slip 
bands which nucleate at the crack tip in the LJ simulation are clear signs that the dynamics 
of the plastic response and the resulting propagating shear modes are crucial aspects of the 
problem. 

Figure |5| shows one example of a local region before and after rearrangement. This 
rearrangement took place in the early stages of the ductile simulation prior to significant 
blunting a small distance in the y-direction from the tip. The arrows denote the sense of the 
externally applied shear in this region calculated by knowing the asymptotic stress field near 
a crack tip. The figure illustrates that these regions appear to be of the type discussed by 
Spaepen as "flow defects''^ or in other contexts as "shear transformation zones" . That is, 
the region seems to consist of roughly 10-20 particles, the rearrangements seem to be local, 
and the "defect" is not mobile in the same sense as a dislocation. Srolovitz, Maeda, Vitek 
and Egami established that these rearranging regions correspond structurally to "r-defects," 
regions of anomalously high local shear stress.!'! The following section will further explore 
why these regions are "r-defects" and how both the high local stresses and deformation 
dynamics arise from the particulars of the intermolecular potentials. 

IV. ANALYSIS 

These simulations beg the question as to why this change in interatomic potential leads to 
the observed change in ductility. Unfortunately, the current lack of a detailed understanding 
of the microscopies of plasticity in non- crystalline materials will rule out a detailed first- 
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principles theory of such transitions at this stage. However, some important connections 
can be made between these simulations and current theories of dynamic fracture and the 
theory of viscoplasticity in amorphous solids presented in FL.@ In addition, toward the 
end of this section a simplified model of the molecular rearrangements at the heart of the 
viscoplasticity theory will be detailed. This model serves to illustrate how a first-principles 
theory may eventually be developed. 



A. Macroscopic: brittle-ductile behavior 

We begin by considering the theory of high strain-rate crack growth proposed by Freund 
and Hutchinson.0 In this theory the plastic strain rate is considered negligible below some 
shear stress <Jfi ow and above this stress the strain rate ef rises linearly. 

if = e t + e ((J s - a f i ow )//j (4.1) 

Here /i is the shear modulus, <j s is the applied shear stress, e t is the flow rate at yield and Eq 
characterizes the strain rate sensitivity. Furthermore e t <C io and the effect of e t will not be 
important for the purpose of this analysis. Using an assumption of strain rate dominance, 
the theory finds that the energy release rate of the crack is velocity dependent. Furthermore, 
the energy release rate of the brittle crack diverges at both high and low velocity. Between 
these two diverging limits there exists a velocity at which the energy release rate of the 
crack is a minimum. According to this model the crack cannot propagate when driven at 
less than this minimum energy release rate. The value of the minimum energy release rate 



depends on the specifics of the plastic response described in Eq. ( [4.1| ) and, by Freund and 
Hutchinson's analysis 

^«1+C^-, (4.2) 

Wip 'flow 

where G c tip is the bare fracture toughness near the tip, and C is a proportionality constant 
which depends on the shear modulus, density and G c ti . [NB: We will ignore a second term 
proportional to i t /io for reasons discussed above.] 
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In the context of this theory we can ask what would cause one material to propagate a 
brittle crack while another admits a more ductile mode of failure. Since a given mode of fail- 
ure can only result if a propagating solution exists, we can conjecture that the ductile failure 
mode results when the propagating brittle solution becomes, for some reason, inaccessible. 
For brittle behavior to have resulted from the narrowing of the inter-particle potential then, 
the minimum energy release rate should have decreased when the potential well width was 
narrowed. This further implies that the narrowing of the potential either caused a decrease 
in io, the sensitivity of the strain rate to a change in applied stress, or an increase in <Jfi ow , 
the critical stress for appreciable plastic flow. 

The simplest explaination for the increase in brittleness would be that the narrowing 
of the potential raised the critical stress for plastic flow. This is not the case. Bulk mea- 
surements of a fi ow obtained by simulating the two systems in periodic boundary conditions 
with zero applied pressure and a constant applied shear strain rate reveal no significant 
differences. &fi ow ~ 0A(e/dl) for both systems. This implies that a change in the flow stress 
is not the cause of the embrittlement in the simulations presented here. 

Returning for a moment to the Freund and Hutchinson model, we note that having 
eliminated <Jfi ow as the responsible parameter for the change in the mode of failure, we must 
consider the parameter eq. This parameter corresponds to the sensitivity of the strain rate 
to an applied stress above the flow stress, essentially an inverse viscosity 

de p J 



e = fi 



da., 



(4.3) 



flow 



In order to explore why such a change in io might arise we will now consider a somewhat 
simplified version of the theory of viscoplasticity in amorphous solids developed in FL. 



B. Mesoscopic: viscoplasticity in amorphous solids 

In the model of viscoplasticity discussed in FL the plastic flow is both rate and history 
dependent. The history dependence of the model enters through a set of state variables 
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which describe the density of "shear transformation zones" (STZ) of the type described in 
section [Til BL These STZ's are theorized to be essentially two-state systems and are assumed 



to have a definite orientation. That is to say that STZ's that are particularly susceptible to 
deformation under one sense of shear may not be susceptible to another, and when an STZ 
undergoes a transition it changes orientation so as to be susceptible to an opposite applied 
shear stress. For the sake of simplification the STZ's are assumed to be either perfectly 
aligned with the applied stress or anti-aligned with the applied stress. Furthermore, the first 
few non-linear terms which describe the dynamics of these STZ's are conjectured. These 
terms represent an assumption that inelastic work done on the system may generate new 
regions or eliminate existing regions. 

In FL the rate of plastic strain is related to the rate at which STZ's transform between 
their two states, 

if = V z As [R+ n + - R_ n_] , (4.4) 

where V z is the typical volume of a region, Ae is the increment of local strain due to an 
individual transformation, n± are the population densities of STZ's in each of the two states, 
and R± are functions of the stress describing the rate at which transitions occur between 
the two states. The evolution equations of n± are written in terms of a master equationiH! 

n± = R T n T - R ± n± - C x (a s if) n± + C 2 (cr 8 if). (4.5) 

where C\ and C2 are constants associated with the non-linear terms which determine the 
rate of STZ annihilation and creation. The equations of motion have two steady states: a 
"jammed" or "hardened" state below the critical stress for plastic flow and a flowing state 
above this stress. For the flowing steady state 

n± = % * kaFcW (4 ' 6) 

The specifics of the choice of the functions R± and their dependence on the stress are 
important for determining the time dependence of the plastic flow. In FL the transition 
rates are written as volume activated processes. That is, the rates are written in the form 
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R± = Ro exp 



AV*(±a s 



(4.7) 



where for the purpose of this analysis we will assume Ro to be a constant attempt frequency, 
Vf is a free volume per particle, and AV* is a free volume needed to activate a transition. 
The volume needed to activate the transition AV™ is a function of the applied shear stress 
which is chosen to have the simplest one parameter functional form for which the volume is 
assured to be non-negative. 

AV*(a s ) = V* exp(-(7 s //2) (4.8) 

where V Q * is the free volume needed to activate a transition at zero stress, and Jx is a modulus 
which characterizes the sensitivity of the activation volume to the applied stress. Note that 
in general V * 3> Vf and these rates are negligible unless a s ~ +/L Since we are considering 
the material response around a s = +afi ow , we are in a regime where R + ^> 

Taking this formulation of the transition rates into account, we can consider the rate of 
deformation described by Eq. ( |4.4| ) in the the steady-state flow regime of Eq. ( |4.6|) . 

if « V z AsR + n + = ^-[aj^ - a; 1 } (4.9) 

where aji ow = (C% V z Ae)^ 1 . So, we can evaluate Eq in Eq. (|4.3| ) using Eqs. (|4.7| - |4.9| ) to be 



£ o = TT^i — exp 



V* 



(4.10) 



v .f 

This last equation provides a first clue as to which aspect of the microscopic behavior is 
responsible for the observed change from ductile to brittle failure. First, note that we can 
reasonably neglect the prefactors to the exponential since the effect of changes in these terms 
will be substantially less dramatic. Furthermore, the ratio V^/vf, which was already noted 
to be a large number, is expected to depend primarily on the relative sizes of the particles 
which are the same in both systems. Since the possibility of a substantial change in <Jfi ow has 
been ruled out, the only remaining parameter in this expression is \x. The double exponential 
causes e to be suppressed by a factor of exp(— V */vf) when JJl becomes large. Moreover 
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Eq. (|4.10|) is most sensitive to changes in Jx when /i « Ofi ow . The investigations of analogous 
amorphous systems via computer simulation in FL suggested that JJl does indeed fall in this 
range. 

The observation that Eq is highly sensitive to changes in ft means that the sensitivity 
of the material flow rate to a change in applied shear stress is highly dependent on the 
sensitivity of the deformable regions (STZ's) in the solid. Relating this to Freund and 
Hutchinson's fracture model this implies that the observed change from ductile to brittle 
failure seems to be due to a corresponding change from "floppier" to "stiffer" weak regions 
in the solid. But while the viscoplasticity theory leads us to these conclusions it does not 
elucidate how one might quantify these ideas and relate them to the molecular potentials. 

C. Microscopic: simplified model of a two-state region 

This idea of "floppy" or "stiff" STZ's can be made more meaningful by considering a 
simplified model of molecular rearrangements. The model should be consistent with the 
observations of Argon, Spaepen and co-workers, i.e. it should be capable of rearranging 
in a local way,iHil and also with the observations made by Srolovitz and co-workers, i.e. 
it should be a "r-defect."H@ Figure |6] shows the most stripped-down model of what these 
two-state systems must look like on the molecular level, four molecules interacting via a two- 
body interatomic potential. Since this unit is embedded in the solid it is constrained from 
undergoing translation or rotation. For particular choices of the interatomic potential this 
four-molecule unit is inherently a two-state system. That is to say that for a Lennard- Jones 
or similar potential the energy is minimized by having as many bonds near the equilibrium 
bond length as possible. In this system there are two degenerate ground states, illustrated 
in Figs. |^(a) and §|(c), in which five of the six bonds are of this length. 

Because transitions between the two states of our four-body unit are associated with the 
development of strain in the solid, the material response must depend upon the rates at 
which these transitions occur. At high temperatures the transition rates are dominated by 
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rare thermal events which occur only as exp(—AU/kT), where AU is the energy difference 
between the ground state and the saddle point illustrated in Figure §](b). This is exactly 
the approach used to describe the time dependent strain in theories of deformation kinetics 
such as those of Eyring, Spaepen and Argon. EISiEl Such a formalism is not of use here. 
In the present system AU is on order unity while kT is three orders of magnitude smaller. 
Physically this means transitions will be driven rather than thermally activated. This is, 
of course, cause for alarm. The statistical approach of the theories of deformation kinetics 
in high temperature systems utilized the statistical nature of the energetic fluctuations to 
discern a time scale. How can a statistical theory be developed when these fluctuations are 
not relevant? Instead, consider the solid to be composed of an ensemble of these two-state 
systems, some small fraction of which are close to a free volume induced transition. By 
using this ensemble picture it is possible to preserve the rare-event aspects of transition 
state theory in order to extract a relevant time scale. 

The concept of a free volume induced transition has been mentioned here, but although 
this concept has been discussed in some detail in FL, the molecular details of such a transition 
have not been discussed. In particular, it is necessary to be able to calculate AV*(a s ), 
the shear stress dependent free volume needed to activate the transition. The following 
paragraphs will attempt to answer the following question: For the model two-state system, 
constrained by its surroundings to a certain area, what is the maximum shear that it can 
support before being driven into a different state? If this question can be answered, then, 
given some applied shear stress, it will be possible to determine the free area (our two- 
dimensional equivalent of free volume) at which a region will become unstable. 

With this picture in mind let us consider in some detail what is going on physically. We 
can parameterize the energy of the four-particle system by only two parameters, its x and y 
dimensions, if we constrain it from rotating, translating or deforming in an asymmetric way. 



Here U can be any two-particle potential, but we will concern ourselves with Ulj and Uclj 




(4.11) 
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described in section [II A| . Furthermore, consider the case when the area (two-dimensional 



volume) of the system remains constant by imposing the constraint A = xy. We can define 
a local equivalent shear stress 

Z s (x,y) = -(-- — . (4.12) 

2 y ox x ay 

At this point it is possible to understand why such two-state STZ's would be visible as 
"r-defects," the regions of anamolously high shear stress described by Srolovitz, et. a/.i'@ 
Consider the condition for the lowest energy of the configuration, 

dU A ldU ldU fAin . 

^ = + (413) 

where the path of constraint is traversed by a unit speed curve parameterized by s such 
that ds 2 = dx 2 + dy 2 . We immediately note that the condition for equilibrium is not 
the same as the condition for zero shear stress. In general these two conditions are not 
simultaneously satisfiable. It is important to note that an exception to this, i.e. a case in 
which the lowest energy configuration has no shear stress, is the case where the molecules 
interact only via nearest neighbor interactions. This is particularly interesting in light of the 
simulations since the limit where A 3> 1 in Eq. ( |2.3| ) is the limit of solely nearest neighbor 
interactions. Therefore, we expect that the CLJ potential (A = 1.5), which resulted in 
increased brittleness in the simulations, should show lower levels of internal shear stresses 
and fewer "r-defects" than the LJ potential (A = 1) which resulted in increased ductility. 
Thus, the microscopic model strongly suggests that the range of the intermolecular potential 
is crucial in determining whether these STZ's are visible as "r-defects.". 

We return now to the question of when the two-state STZ will become unstable to an 
externally applied shear stress. The condition for instability can be written 

ds + y ® x vy 

where we have again traversed the path of constraint by a unit speed curve. We can now 
define the equivalent of a free volume in our system, Af = A — Aq. Here Aq is the equilibrium 
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area of the four-body system at zero applied shear stress ~ \/3 (g?q). Figure [7| shows the 
value of Af(T< s )/Af(0) at which instability sets in for values of the shear stress S s . In order 
to illustrate the suppression of AV*, and by analogy Af, as the stress is raised to the flow 
stress, Figure ^ spans a range of shear stresses up to our observed flow stresses in the fracture 
samples. This graph looks suggestively similar to the form guessed in Eq. fl4.8Q . In actual 
fact, however, in the vicinity of S s = 0, Af is a power law and not an exponential decay. This 
suppression of the activation volume can be related to values for p, in FL. The activation 
area at S s = 0.4 is 58% of its value at zero in the LJ case and 74% in the CLJ case; this 
corresponds to pij = 0.73 and Pclj = 1.3. Thus, longer range intermolecular potentials 
correspond to a solid with "floppier" two-state regions. This result is significant since longer 
range potentials also implied larger local stresses, and the existence of "r-defects." As 
expected from the previous analysis, the toy model with the CLJ potential has a higher 
value of p, and, therefore, corresponds to a solid with "stiffer" two-state regions. This is in 
keeping with our expectations since "stiffer" two-state regions should also correspond to a 
lower value of Eq and, therefore, by Freund and Hutchinson's model to a lower minimum 
energy release rate for brittle fracture from Eq. fl4.2|) . The CLJ solid is observed to undergo 
brittle fracture. 

The analysis presented here is clearly only a first step toward a rigorous first-principles 
theory of brittle-ductile transitions in noncrystalline solids. Future investigations will hope- 
fully allow more explicit connections to be made between the molecular level structures 
which are quantifiable via diagnostics such as D^ in and the observed fracture behavior. 
Progress will require further developments in our understanding of the molecular physics of 
deformation in non- crystalline solids. 
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FIGURES 




Interparticle Distance 

FIG. 1. The LJ and CLJ potentials. Energy is given in units of e. Interparticle distance is 
given in units of do • 
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FIG. 2. Stress averaged throughout the sample versus time for the CLJ (solid) and LJ (dashed) 
simulations. The higher stress for the onset of fracture in the LJ case implies increased dissipation. 
Stresses are given in units of the critical stress for failure of an ideally brittle material with the 



same elastic properties. Time is given in units of doy/m/e. 
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t = 400 

FIG. 3. Frames from the CLJ (left) and LJ (right) fracture simulations. In each set the frames 
are shaded by the parameter -D^m defined in Eq. (|3.1|). Dark regions have undergone the highest 
amount of non-affine rearrangement. The shading saturates when -D^j n = 1. 
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t = 600 




t = 800 




t = 1000 

FIG. 4. Frames from the LJ simulation showing the nucleation and growth of a void in the 
vicinity of the crack tip. The frames on the left are close-ups of the crack tip. The frames on the 
right are shaded as in Figure ||. 



25 




t = 40 t = 50 

FIG. 5. A local region before and after non-affine rearrangement. The molecules are shaded by 
Dmim * ne a moun t of non-affine rearrangement. The arrows denote the approximate direction of 
the externally applied shear. The ovals are included solely as guides for the eye. 





(a) (b) (c) 

FIG. 6. Diagram of four identical interacting particles making a transition from one stable 
configuration to another. The middle configuration is the saddle-point configuration for this tran- 
sition. 
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FIG. 7. The activation free volume (area) at shear stress T, s divided by the activation free 
volume (area) at zero applied shear stress, Af(Ti s )/Af(0). The activation free volume (area) 
corresponds to the excess volume at which the four particle system illustrated in Figure [] becomes 
unstable to the applied shear. Stresses are given in units of e/dp. 
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